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ABSTRACT 

The use of radar tracking in aircraft intercept operations intro- 
duces a random noise factor into the control system. Some techniques 
have been proposed for organizing and smoothing the radar data prior to 
insertion in the control system, however the efficiency of these tech- 
niques in providing an optimum interpretation of radar data is suspect 
in the case of maneuvering targets. A simulator design using a general 
purpose digital computer is proposed to enable observation of the behav- 
ior of certain track smoothing methods in a typical intercept problem. 
The method of aircraft vectoring control in the simulator design is dis- 
cussed, The effects of angular and linear target acceleration on some 
tracking methods is described and illustrated. Application of sampled 
data theory to digital computer programming is shown in a simple filter 
design. Simulator was programmed by author and the USNPGS CDC 1604 com- 


puter was used to conduct testing. 
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1. Introduction 

The use of radar as a data gathering device has gradually been sup- 
plemented by the introduction of processing and computing equipment de- 
signed to make the most efficient interpretation of the information 
thusly obtained. One factor that complicates the process of interpreta- 
tion is the presence of random error in the data which results from the 
resolution qualities of the radar. The problem of filtering this data 
has been the subject of many papers since that presented by Wiener two 
decades ago. [1] In some applications, however, the theory presented 
cannot be directly used. One such situation is that which remains a 
primary problem in the air intercept control system; the tracking of a 
maneuvering aircraft. Several methods of tracking and smoothing have 
been advanced for use in this system which purportedly provide suffi- 
ciently smoothed data for intercept geometry calculations. Whether or 
not these methods, using weighting function response techniques over a 
finite time, can provide data of sufficient smoothness and accuracy is 
a question which can be answered only by controlled testing in a typical 
system. The most practical way of conducting such tests is through use 
of a simulator which will duplicate the operation of the entire air in- 
tercept control system. 

The purpose of this paper is to report on the design and operation 
of such a simulator, programmed into a general purpose digital computer, 
and to investigate the results of some limited tests conducted thereon. 
The programming methods and subroutines used in constructing the simu- 


lator were chosen with little regard for arithmetic efficiency or mini- 








mizing storage space, and hence are not submitted as system proposals. 
Rather they are presented as being representative of methods which might 


be used in such a system, 








2. Preliminary Analysis 

The data obtained from a radar set provides range, azimuth and ele- 
vation signals in discrete, sampled form, with sampling rate of 1/T 
where T is the period of antenna rotation. A video processor operating 
on this output would ideally remove all spurious data and pass on target 
information with a blip to scan ratio of one, still in discrete form. 
For simplicity it will be assumed that the miss probability is zero. 
This sampled data, which includes the target position signal plus a ran- 
dom error function in range and azimuth due to a position uncertainty 
inherent in the radar, must then be operated on by a position and velocity 
smoothing procedure to obtain a best estimate of track of each target. 
Once the tracks are established, computations are made based on desired 
tactics so as to control the interceptor aircraft into a particular 
position relative to the target. The maneuvering of the interceptor in 
response to the control commands is interpreted by the radar and the 


system is then closed except for the inclusion of human decision factors. 
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Figure 1. Block diagram of Simulator System 








The simulation of this closed system is accomplished using a digi- 
tal computer (CDC 1604) with FORTRAN chosen as the programming language. 
The program runs in fast time while maintaining a real time index for 
scaling velocities and turn rates. Of the two possible methods of con- 
trol computations recognized, that of performing a complete solution at 
the onset and updating as necessary, and that of making a partial solu- 
tion with each sampling, the latter was chosen to permit overall use of 
sampled data theory without the need of adjusting for a variable time 
base. A flow diagram of the simulator design used is included as Ap- 


pendix I. 


A. Intercept Control Computer Program. 

The method by which the computer generates control information 
is not considered to be a emtecall factor in the simulation of the system. 
However, it is felt that a general knowledge of the method used will as- 
sist the reader in following developments made later. 

The geometry of the intercept was chosen to be one in which the 
interceptor is vectored for a collision course intercept with an offset 
point which is a function of target crossing angle, velocities of target 
and interceptor, and a particular bank angle to be used for the turn 
from the offset point to the target heading. This offset point is cal-~ 
culated so as to put the interceptor at a point one half mile astern of 


the carget on completion of the final turn. 
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Figure 2. Offset Point Geometry 


In Figure 2, the X and Y relative distances are computed using the 


formulas; | 
_ Ge | 
Kee, = Ve ROT. .& - RAD, Sin (1621) 
Yee. = RADz (1-005 (/Gp!)) 
, 
as | Gel 


and, by flagging the side from which the intercept is to be made, these 


relative distances are rotated to the x, y reference by 


Xaes = Xeer Cos G. - SIDE: Veer Sin G, 
Yass = kee, sin Op + SIDE: Voes Cos Gy 
During the phase of the intercept where this point is being intercepted, 


the offset coordinates are added to the target coordinates to generate 


a pseudo-target on which the collision course calculations are made. 





The calculations for the collision course intercept are based on 


relationships in Figure 3, 





Figure 3. Intercept Geometry 


where the bearing to target is SRG = tar" (Ye Ye Vike - ke\ | 


and the relationships 
GP fh a = 77 + BG 


are used in the sine law to find HDG solution 
[ sel 
and intercept heading then becomes H40G= BRG+). 


No mention of wind is made in these calculations. Since wind at 
normal interception altitude may be as high as 25% of the aircraft spe- 
eds, considerable difference will exist between actual track fed in by 
the video processor and computer output commands. Also, target heading 
and true air speed will not be coincident with track parameters and will 
cause errors in determining a precise intercept course and offset point. 


Although it would be impossible to predict winds throughout the inter- 








cept progress for complete compensation of these effects, a good esti- 
mate could be used based on prior knowledge which would reduce these 
effects to an insignificant size. In the program used, there is a 
routine inserted to convert the track smoothing output into a no-wind 
situation prior to insertion into the computer program, thereby elimina- 
ting the necessity of including wind in the actual intercept computations. 

The intercept, on reaching the offset point (within distance limits), 
then proceeds with a constant rate turn of ra) = 30° toward the target 
heading. To realistically complete this turn, when the relative target 
bearing reaches a computed antenna train angle, the turn changes into a 
target tracking course which simulates the conversion to aim dot ste- 
ering from the airborne computer system. 

The output of the intercept control portion can either be continu- 
ous, as in the case where data link communications exist between computer 
and interceptor, or of the step output form, where commands are passed 
only on significant change of the computer solution, which duplicates 
the present manual control method. The former method involves the gen- 
erating of a command acceleration (linear and angular), and the latter 
the setting of a heading and altitude. The latter is used in this pro- 
gram, and the method selected for generating the command headings and 


altitudes will be investigated later. 


B. Interceptor Response Simulation 
The reaction of the simulated interceptor to commands from the 
control computer section is simultaneous. Although pilot reaction time is 


not negligible, and would in practice be accounted for, it is assumed 








to be compensated for in this simulator. Hence, on command for a turn 
to a new vector, the simulator sets the turn rate corresponding to the 


present true airspeed using the aerodynamics relationships, 


RAD = 189.3 (TAS) / tan 


ROT =(0.00529 tan 6)/TAS 


where PHI = _ bank angle in radians 
TAS = true airspeed in miles/sec 
ROT =~ rate of turn in radians/sec 
RAD = _ radius of turn in miles 


The turn rate is then converted to radians/time base, and the x,y co- 
ordinates of the interceptor are modified according to the relationships 


shown in Figure 4. 
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Figure 4. Turn Compensation 


where (epee G, +wt ; (t= radar scan time ) 
where the arc length is the actual distance travelled, 


S= RAD (wt) 
and the chord length, 
a ult 
d= 2RAD sin(4 ) 


are related by the factor 





we 
sin(>) 2 
Cony eee = é = L a rere eae 
s wot / om PR Han itayaee 
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and advance of the coordinates are made by the formulas 


Mee, = Xp, +|4 (Conv) Cos (Sc + 2) 4 We] (i 

Ye nes = Yen + le (Conv) Sin (Ge + ot) aw] fe 
where WX and WY are wind drift inc ements. The target is advanced in 
the same manner according to a programmed maneuver. Headings are ad- 


vanced according to turn rate, 


ay aes Ge, + a,t 
Gr nay : Gr, +W,t 


and when the heading comes within a present range of the mew vector, a 
random deviate is looked up and the interceptor course is set within 
+ 5° of this vector (normal distribution) to simulate error in heading 
control. 

Speed and altitude control responses are simply straight line vari- 
ations, with speed change set with a constant acceleration of + .00028 
mid @ey/sec™ » and altitude change based on the need for reaching the new 


altitude on expiration of 50% of the time to go until intercept. 


- 


C. Radar and Video Processor 
The ‘actual’ blip positions representing true target positions 
are taken from the interceptor response simulator in rectangular coordi- 
nates and converted to polar form. According to the resolution capabili- 
ties of the radar set being simulated, a Gaussian noise is added to the 
range and bearing signals, which are then translated back into rectangu- 
lar form. The processor output is set at a blip to scan ratio of one by 


this procedure, and the hit probability is assumed to be one. 








3. Systems Design 

A. Control Methods 

The intercept computer generates a solution based on present 

geometry with each iteration. Since it is assumed that this simulator 
is to represent a system which is limited by the computer-interceptor 
communications link to a step command output, a non-linear system must 
be considered in generating this output. The simplest system is ome 
wherein a new heading command is given at the time the computer solution 
differs by a set amount from the last command, the new heading command 
taking the value of the solution at that time. If the interceptor were 
able to respond instantly to this command, and hold the exact heading 
ordered, with no time variable parameters involved, this ome order would 
be correct during that particular phase of the intercept. Such is not 
the case, however, since the response time of the interceptor is pro- 
portional to heading change and turn rate, and since parameters are in- 
volved which are time variable, such as linear acceleration to intercept 
speed and target heading during track smoothing. Also, there is am in- 
herent steering error and speed control error on the part of the inter- 
ceptor which is not predictable and which along with erroneous wind 
estimates causes a drift in the computer solution. Hence, by the time 
the interceptor reaches the new heading, it is mo longer the correct 
heading, and the result is a sucession of heading commands lagging the 
correct value. 

One of the operational factors involved in controlling the inter- 
ceptor is the necessity of reducing the number of control commands to a 


minimum, to avoid communications saturation. To accomplish this end, it 
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is necessary to anticipate interceptor response and computer solution 
behavior in determining the next command to be given, in order that a 
more nearly correct heading will result (the event in which the inter- 
ceptor will take the exact heading required has a probability of zero). 
A simple method of prediction is that of using the derivative (continu- 
ous case) of the required heading and assuming a linear change during 


the time the interceptor completes the rurn, as in Figure 5. 






-——-WNEW VECTOR ASSIGNMENT 
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/ heading 


HEADING 


TIME 


Figure 5. Sketch of Derivative Prediction 


The error involved in this method depends on the amount of non- 
linearity present, which is indicated by the order of the polynomial re- 
quired to fit a curve to sampled points. Using data from an early trial 
run it was found that a good approximation could be obtained by a third 
order polynomial at short ranges, and a second order polynomial where 
the range to target was over 50 miles. This would indicate that the 
simple constant plus derivative (difference) method might deteriorate 
at short range where heading control is critical. 

Another possible predictor, of greater complexity, would be an ex- 
trapolator using a difference table of some chosen order. Im order to 
compare this method with the previous to determine if any signifi. ant 
advantage might be obtained that would outweigh the computer time in- 


volved, both methods were synthesized in the simulator. The main ad- 
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vantage in finding an efficient and accurate predictor lies in the re- 
sultant minimization of steering error afforded by providing the most 
correct mean heading about which the steering error is applied. 

To predict the required HDG at some future time based on past his- 


tory, a third order polynomial is fitted to the last four HDG solutions, 


using the divided difference table. 


a, 
to f, 
fh 
f Mf, 
Af : 
te io } Nt A £, 
Af. | 
ts fs 


Although the actual function is of higher order, the error in assuming 


a third order fit is small compared to steering errors. Hence it will 
be assumed that 
Ar4, = A? f, 
Now, rather than using Newton's forward interpolation formula, an extra- 
polation is made based on last differences, 
5 oe 3 3 
Af, = A4,+Af, i Af; 2 Af, + af, 
In order to make this system adaptable to varying time base, the special 
formulas are used, 
2 2 3 zs 


and to find f, at some given time following tz 9 
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ie . ie in (ty-ts)} Aft (ty-te)] AF + (ty-t,)o°%4o]} 


where oh is solved for using 


ae d(f@)) acer ae ( f5 — Peesenr HEavwe) 
4 
atte: (to- Ate) 
according to the sketch shown in Figure 6. 
f- — — — NEW VECTOR 
fi ASSIENMENT 


HEA OING 





interceptor turn 
rate > ROT, 


REQUIRED 






Figure 6. Sketch of Extrapolation using Last Differences 


On initial vector assignment, extrapolation may cover a wery long 
period and involve a significant error, however, since the solution is 


made continuously throughout the resultant turn, it is self-correcting. 


B. Tracking and Smoothing Procedures 


The output of the video processor consists of sampled signal 
plus noise information, in the sense that noise exists as a random 
error in position, with predictable variance and an assumed Gaussian 


distribution. In the translation of these polar form data into rectan- 


gular coordinates, the independent random variables associated with 


range and azimuth error become functions of position. Assuming that 


range and azimuth signals have the respective deviations Se and % 
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then in converting to the rectangular form, the variance associated with 


the x and y coordinates become 


is Z Pipe ost? 
QF G tos @ + OP Sin O 
and 
7 og O 2 2 2 
OG, = Gp Sin® + Op cos © 


Graphically, the joint probability distribution of the random variables 
in polar coordinates at a particular point would appear as in the figure 


below, 





Figure 7. Sketch of Joint Prob. Distribution of Typical Radar Report 


with the transformed x-axis random variable distribution shown by the 
cross-section. On successive scans of the radar, with the processor 
providing a blip to scan ratio of one, a set of random variables is pro- 
duced whose correlation is dependent on the shepe of the actual track. 
Given a sample of sufficient size, a regression curve could be fit based 
on measurements obtained using statistical analysis. Im the event that 
the dependence of x and y were linear, as in a straight track, the slope 
and y-intercept could be found directly. However, this method is not 
satisfactory since the information required is available only after 
future samples are taken, and since new position information is not di- 


rectly available. A more desirable means of estimating the linear re- 
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lationship of the x and y coordinates would be a recurrent updating 
of a track estimate based on previous values using a weighing function 


according to track length. 


(1) Constant Parameter Track Smoothing 
Some of the early work done in the field of track position 
and velocity smoothing is that presented by Boxer in extention of the 
formulas of S. Thaler [2]. He usis the difference equation representation 
of a digital computer operation to construct and analyze a feedback net- 
work representing a positional tracker with velocity aiding. In the 


block diagram, 





Figure 8. Block Diagram of Velocity Aided Positional Tracker 


| 
the Za term represents a digital accumulator and the z in the numer- 
ator of the other term provides output without delay. The overall trans- 
fer characteristics of this system can be determined by z transform 


methods, yielding in recurrence form: 


Yo (4427) + dvs (teT)- BVO(4)= (tps KAGE) + Ky GY 
where 


. 
oe - ¥ 
Aé)= @(t+7) -V) Ae Xpo! Q= Kot Ky mee Te 4 
T 
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x 
Letting the driving function v; () be equal to zero provides the char- 


acteristic equation, whose Laplace transform is 
x ‘ ae 
Re (2) B+ X2-BT = 0 
2 
The stability criterion here is that & and 8B , the roots of #4 +z 


im 6 = O , have absolute values less than unity. This describes the 


allowable values for Kp and Ky according to the figure 


4 


ky 


e) > Kp 
Figure 9. Stability Area for Smoothing Parameters 
to lie within the shaded portion for stability. Further analysis of 
stability is advanced by Boxer in forming the over-all transfer function 
of the system and taking the denominator as an equivalent Nyquist equation. 


The same results obtain. 


By labeling the signals as shown in the figure below, 


a see ae 





Figure 10. Block Diagram of Velocity Aided Tracker 


we can develop the relationships as follows: 
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0 
Ae aye Neen, 


from which we get 


/ 


0,@)= K2@)+4Oe@2 
or, 
Oft)= WM L@)+ 0-7) 
Xe = ky ele 
where the input is taken to be (X,-x,) . The input to the second ac- 
cumulator is then 
ie he EP 


and the transfer function 


«(e) (2) 5 Qb)=2,4-7)+ 04-7) 


(ee 
gives the output as 
A ~—. ivy aA a 
Xn = Kner + kp Ce) ~ Xne) 


If we define the quantity 


oe A ~ a 
Xne a Xnes t Kp (vas -Xp-1) 
then we can write 


ie A rs 
Xn = Xpy-, t Kn- 


Thus showing that this development leads directly to the prediction and 


update formulas generally used for track position and velocity smoothing. 


- A ho aA 

Xn = Xn +A(Xp- Xn) (1) 
oe = om +8 (Kn - kn) cee for 7= 19 (2) 
Xno1 = Je ee 
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where the substitutions X= Kp and B=k, have been made. 

The preceeding plot showing the stability limits for & and a 
indicates that there is no unique set of values, and the possibility is 
that the performance of the system may be adjusted for optimum response 


to a particular set of input values by varying these parameters. 


(2) Variable Parameter 
By defining a variance in the input signal, and using the 

least squares method of fitting a straight line to the input, FCPC has 
developed recurrence formulas to be used in a computer to control the 
Qa , 6 parameters [3]. Such a method, when used with an appropriate 
turn logic, provides optimal smoothing during those portions of the 
tracking problem where no turns occur. The derivation of these formulas 
follows from the definition that the best value of a set cof measurements 
is that for which the sum of the variances of the observed values is a 
minimum. Thus from (1), 

- ~ A 

Kn = Xy Xn t (1-Mn)X, 
we define the variance of x, to be the square of the constants of the 


A a 
independent terms Y¥, and X, times their respective variances, 


Vn = x, Va a (/- %) Vy 
and from (2), 


~ 
Kn = Bn Kn = Bn X,, + hig 


using same definition for variance of x, 2 


= % cf cs 
Vin = Bn Vn Gr An Viney + (i- Bn) Vner - 2 Bp (/- Bn) COV,_, 


where covariance COV is the expected value of the product of the deviations 
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of the variables from their expected values. 

The FCPC report develops a set of final equations based on this type 
of analysis in which the values of & and & are computed each pass 
according to equations that would minimize Va and Vn respectively for 
a constant velocity input variable [3]. These equations are not repeat- 
ed here due to the classification of the reference, however, they are 


used in the system and identified as the variable parameter method. 


C. Data Filtering 
Once the program determines smoothed valves for target position 

and component velocities, X 6 Y 5 me yy. the track is calculated and 
ground speed computed directly from ¥ and y . These are then fed in- 
to the computer for intercept calculations. It was noted that small er- 
rors in heading and track has a multiplicative effect in the offset 
point calculations, causing significant excursions from the expected 
values. By comparing the offset point track in Figures 13 (noiseless) 
and 16, it can be seen that this effect is detrimental, causing exces- 
sive hunt for heading on the part of the interceptor. It was thought 
that a special smoothing technique in the computing program itself, in 
the form of a sub-routine, would be beneficial to stabilize the calcula- 
tions. 

Using the relationships which exist between the sampled data system 
and the digital computer, and with the program arranged to iterate each 
radar sweep, it is possible to apply sampled data system theory directly 
to the problem. Hence, we can construct a low pass filter, using the 
analogy of a friction damped mass in a simple position feedback network. 


Such a system can be drawn as in Figure 11 as a continuous system. 
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Figure 11. Simple First Order Positioning System 


Its sampled data equivalent is shown in Figure 12. 





Figure 12. Sampled Data Equivalent of Figure 11 


where the @ and XK can be chosen to obtain the desired ‘filter’ char- 
acteristics. 
Since the noise in the system is predominantly of the frequency 
a » the cut-off frequency of the filter should be well below this 


value. Writing the continuous signal over-all transfer function 


K 
s*4as7k 





Ro 
WG) = > (5) = 
¢ A © 
and arbitrarily choosing a damping ratio of about ,7 and a baidwidth 


of about Z radians/sec. 
$T 


Ke uw)” = (ZY = (0.196)2 0.0375 (r=4) 


and, 


Q= 2? w,, = 0.275 
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with 
AT = Il 3 &= 0.196 
Since overall transfer contains no finite zeroes, bandwidth is given by 


relationship, 
eS 
BW. = Un[/-29%+-Ve-4gtr 4 ge |* 


which in this case is actually w, : Now switching from Laplace to 
z-transform, the open loop transfer is 


B'L6o] = 2°] 8 25] + a ote 


/-27')(/ -.3333 2.) 


making the overall transfer function 


el _ _ 0.0466 27’ 

R /~ 42867 27 + 0,3333 27? 
Taking the inverse transform of this equation and putting it into the 
time solution, 

C(t) = 0.0466 R (4-7) + 1.2867 C (4-7) -.3333 C (t-27) 
produces a recurrence formula which uses values associated with the two 
previous sampling times in computing output. A delay of TI seconds is 
inherent, and response to various input signal functions is dependent on 
the parameters chosen. It is anticipated that the signal to be filter- 
ed contains no step function components and that the size of steady state 
error to ramp and acceleration inputs is not critical. For the values 
chosen, typical response to a noisy signal with step components is shown 
in Figures 17 and 18. 

The linking of sampled data system analysis with digital programming 


in this manner allows the full use of response, compensation, and stabi- 
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lity theorems in designing the program routine. An important fact to 
consider, however, is that while the system transfer function is arbitrar- 
ily selected, the resultant poles and zeroes are exact and stationary, 
enabling complete cancellation if that form of digital compensation is 
envisioned. Hence, the filtering action of the system would be lost. 
Also, the use of derivative feedback, which merely adjusts the position 
of the control poles, is likewise meaningless since we are free to cho- 
ose at the outset which poles we desire. There remains then the theo- 
rems governing response according to the input anticipated. Since this 
input will likely include step, ramp and aéceleration inputs of some 
magnitude, it seems most desirable to adjust the system for a minimum 
overshoot response to a step input while holding steady state error to a 
ramp input within reasonable limits. Also, it would be possible to set 
initial error at zero with the first step input and to jump the output 
to this value, thus eliminating the initial step response problem. The 
simulator allows us to examine the type input signal generated, and to 
choose the most efficient method of design. The use of the initializing 
routine at the start of the run eliminates initial overshoot problems at- 
tendant with the first step input discontinuity, allowing concentration 
in the filter design of maximum smoothing consistent with allowable 
steady state error to ramp inputs. Redesigning the filter with this ob- 
jective, it is possible to achieve better smoothing without being pena- 
lized by sluggish response to step inputs. It should be noted that the 
object of this filter is primarily to remove random ‘error’ signals from 
system variables in order to stabilize the calculations for the offset 


distances, and that small steady state errors are quite tolerable in 
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view of the magnitude of steering and speed errors already present in the 
computations. Were this not so, then closer attention would have to be 
paid in selecting the velocity and acceleration constants of the filter 
system. 

An alternate filter design using a damping ratio of .9 and 4 =.2 


was used, with step-sense included, having the overall transfer function 


oie Qo8¢g 27! 
R 1 ~/1S17 B71 4+. 2365 2°% 


with a corresponding recurrence formula of 


Ct) = .0848 R(¢-7) + 1.1517 CA-7) -.2305 C-27) 
The logic used for the step-sense, along with the representative recur- 


rence iteration, appear in FORTRAN language as; 


IF (K - 1) 99, 99, 100 


99 Rl = Variable 
Cl = RI 
c2 = Cl 
100 C = Al * Rl + A2 * Cl + AI *® C2 
c2 = Cl 
Cl = C 


Rl = (variable) 


101 CONTINUE 


where Al, A2, and A3 are preset constants. 
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4. Simulator Testing and Analysis 

Three simulator configurations were used in conducting the tests of 
the system. All were programmed with the same geographical positioning 
of the target and interceptor as shown on Figure 13, with headings and 
speeds and altitudes as indicated. In each, the interceptor was program~ 
med to make a beam attack from the nearest side. The program differ in 
that system I has a 'noise' free input, providing exact track information 
to the computer. System II has a constant parameter tracker-smoother, 
and system III has a variable parameter tracker based on least-squares 


analysis. 


A. Control Methods 

Simulator configuration I was intended to allow analysis of 
the control methods outlined in section 3 The results of runs made of 
this configuration are shown in Figures 13 - 15. As the labelling indi- 
cates, in Figures 14 and 15, the simple derivative predictor, when used 
for controlling the interceptor heading, is not able to anticipate the 
rate of heading change that actually exists, and ends up following the 
input at shorter ranges to the extent that an appreciable lag may exist 
at the critical turn point. If the aircraft speeds were increased, the 
input heading drift rate would also increase, causing a marginal situa- 
tion in completing the intercept with the accuracy desired. The extra- 
polator does a better job in making these predictions, and in the pro- 
blem run, resulted in fewer heading command outputs. It should be noted 
that in the testing of these methods, there was a random error factor 


inserted in the interceptor response to each heading command, of the 
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range shown in the figures, so that the results would be more directly 
applicable to the actual problem. In this light, we should also have 

to consider the track accuracy required in the practical problem. Since 
this phase of the intercept brings the interceptor to a point from which 
a final rendezvous type turn is made, during which time compensation can 
be made for a fairly large positioning error, there could be a wide lati- 
tude for error without affecting the final success of the intercept. In 
this analysis, however, the major concern is to remedy knowm sources of 


error so as to keep the cumulative error to a minimum. 


B. Simulated Tracking and Smoothing 

The basic organization of the intercept control computer is 
complex in the switching functions and decision routines that are neces- 
sary to complete a simple intercept. The equations themselves are 
straightforward, however, and the construction of the simulator becomes 
a problem in programming. The operation of this control system when 
linked with a noise-free feedback path of position and velocity as shown 
in Figure 13 is smooth and predictable, due to the linear behavior of 
the variables involved in the intercept computations. Such is not the 
case when the feedback path is impressed with random noise. Now the in- 
put variables to the control unit are non-linear, and this effect is 
magnified in the computer solutions, resulting in the type of intercept 
shown in Figure 16. The random movements of the offset point track 
demonstrates this effect. In this intercept, the fixed parameter method 
of track smoothing was used with a turn logic shown in Appendix II. The 
turn logic is designed to sense a departure of an incremental heading 


from the previous smoothed value by a certain set margin, at which time 
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an initiation of a new track is signalled. It does not sense a departure 
of incremental velocity, however, and this appears to be a serious draw- 
back in using this turn (only) logic. In Figure 17 the response is 
shown of this tracker to a ramp velocity function corresponding to the 
acceleration of the interceptor to attack speed. Since the interceptor 
speed requirement is interpreted from the target speed, which is derived 
in the same fashion, the target maneuver in the vicinity of the 65th 
iteration causes multiplicative effects in the program. By analysis of 
the results shown in Figure 17 it is apparent that this tracking tech- 
nique cannot be used directly throughout a turn maneuver, even by initi- 
alizing with each scan, in view of the error excursions generated each 
time a turn occurs. This is due in part to the error involved in using 
the chord length throughout a portion of the turn as being the actual 
distance travelled. This error is equal to the departure of the Sinus 
curve from unity for the angle of turn during which the approximation is 
applied. A much more significant error is also present, prior to the 
conversion of the x and y velocities into heading and velocity. As 
the turn progresses through the heading of one of the major axes, the 
velocity associated with the other axis changes sign. In the recurrence 
or weighting function, no provision is made to accomodate this rapid a 
change of velocity, and hence the velocity lags while at the same time 
the velocity associated with the other axis remains fairly constant, 
producing a significant error during conversion to resultant heading and 
velocity. On actuation of the turn logic, and initializing of the track- 
ing routine, this error is temporarily resolved. The critical value ap- 


pears to be the limit set in the turn logic which determines when an 
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actual turn has begun. If this limit is too small, input signals with 
large deviations would trigger the logic and upset otherwise good smo- 
othing. If it is too iarge, an excessive amount of turn would occur be- 
fore triggering the logic, resulting in large errors im speed and heading 
during the turn. In the tests made, this variable was set to sense when 
any incremental variation in heading differed by more than 114 degrees 
from the smoothed value, and was not changed during any of the runs made. 
It is also apparent from the figure that the tracker fails to follow a 
linear acceleration, which is also due to the constant and relatively 
small weighing factor B that serves to update the velocity. The de- 
pendence of the tracker routine on the turn logic chosen is evident in 
Figures 21 and 22 where the radar resolution was purposely deteriorated 
to increase the noise in the system. It can be seen that the output of 
the tracker more closely follows target maneuvers in heading and speed, 
at the expense of providing poor smoothing during these maneuvers. This 
resulted from the more rapid triggering of the turn logic by the input 
signals of larger deviation. 

While a re-design of the tracker and turn logic equations is cer- 
tainly indicated, the probability is that the random nature of the out- 
put will not be completely reduced. To investigate the possibility of 
using additional smoothing in the control program, the filter described 
in section 3C was inserted to smooth selected input values. The degree 
of smoothing provided in the case of interceptor velocity is shown in 
Figure 17, and in the improved overall solution as shown in Figure 19, 
where filtering was performed on interceptor velocity, target velocity, 


and target crossing angle. 
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Of interest is the fact that the hypothetical radar selected for 
these tests had relatively superior resolution capabilities of about 1/8 
mile and 15 minutes. A reduction in this performance to a resolution of 
0.3 miles and 35 minutes, using the same scan rate resulted in the solu- 
tion shown in Figure 20, with filter performance shown in Figures 21 and 
a ae 

The third simulator configuration used the tracking and smoothing 
procedures based on least squares fit as outlined in section 3B(2). 
Again, this method is intended for straight line tracking amd required a 
turn logic routine. This routine is shown in Appendix III. In the test- 
ing of this configuration, it appeared that a separate logic was needed 
to sense linear acceleration as well as angular, since the variable 
weighing function associated with velocity would not respond to an in- 
crease in target speed rapidly enough, as seen in Figure 18. Since the 
performance of the tracking methods used seems to be largely dependent 
on the choice of turn logic used, an evaluation of these methods cannot 
be made on a maneuvering target, and no conclusions can be drawn as to 
the merits of these methods in the intercept control system. 

If the tracking systems were considered as being representative of 
some sampled data system, such as the analogy shown in the constant para- 
meter system, with the recurrence formulas representing the weighting 
function obtained by the inverse transform of the transfer function, 
then we could assume the system input to be a noisy ramp input (of posi- 
tion) and a noisy step input of velocity. Initialization of the routines 
corresponds to setting initial error at zero, and the response of the 


system is designed to produce minimum error for the expected input (con- 
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stant velocity). Now if we assume that the position input will be of 
second order, as is the case with a maneuvering target, it is necessary 
to redesign the system for minimum response error to this type signal 
also. Then, realizing the dependence of the x and y coordinates of 
a realistic target, we should also solve for heading and velocity and 
perform identical smoothing on these variables. With such a system, a 
separate turn logic would not be necessary. Since the design of such 
a routine is not the intent of this paper, the foregoing is submitted 
merely as a proposed approach based on the observances made thus far 


using the simulator routine. 
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5. Conclusions 

The use of a general purpose digital computer in simulating an air 
intercept control problem has enabled the author to investigate the be- 
havior of a hypothetical yet typical control system and two proposed 
tracking and smoothing procedures. The results of the limited investi- 
gations made are valuable only in so far as the arbitrarily selected en- 
vironmental conditions are concerned. It was found that the constant 
parameter smoothing method relies heavily on the appropriate choice of 
turn logic in producing a heading output of any instantaneous accuracy 
throughout a turn maneuver, and that velocity output throughout this 
maneuver is quite erratic. The choice of this turn logic appears to hinge 
on the scan rate and the anticipated rate and duration of turn, rendering 
this elementary tracking method inefficient for use in the system. 

The variable parameter method proved to be rapidly convergent to a 
very smooth heading and velocity output during straight line portions of 
the track being followed. However, like the constant parameter method, 
it depended heavily on an appropriate turn logic for even the most ap- 
proximate outputs during a turn maneuver. The basic problem in both 
methods was the inability of the tracker to follow component velocity 
changes throughout a turn or throughout linear acceleration. A turn 
logic which would function in a manner similar to the tracker in weight- 
ing a turn rate variable and linear acceleration variable appears to be 
the next logical extention of these methods, and is suggested as a pos- 
sible area of further study. 

The requirements of the control system in computing a satisfactory 


solution to the intercept problem were found to be quite flexible. The 
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most detrimental effects were seen to originate more from noisy input 
variables than from the erroneous but smoothed, filtered form of these 
variables for the type intercept solution made. If the offset point 
method were to be used with an additional beam displacement, as is most 
probably the case in present air tactics, then the effects of noisy vari- 
ables would be more acute and would accentuate the need for additional 
filtering of the heading and speed variables. 

It is hoped that this paper will provide a partial ground-work in 
the intercept problem from which further investigation may be made of 
some of the problem areas revealed. It is suggested that the program 
shown in Appendix I be oriented about a method of simultaneous display 
and direct computer control in order to allow more rapid and flexible 


testing. 
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APPENDIX I 


BLOCK DIAGRAM OF SIMULATOR PROGRAM 
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MORE NOTES ON SIMULATOR PROGRAM 


1. MODE identifies the type run calleds; 1 for beam type, and 2 


for head-on type. 


2, PUASE indicates-the phasé:of:the intercept; 1 until offset 
point reached, 2 during attack turn until 3 when intercept is 


taken over by aircraft. 


3. STEP is a 20° increment to facilitate faster rotation of the 
bearing line during head-on intercepts, and faster off-track 
displacement during beam runs, It is added to HDG prior to en=- 
try into the heading control and command section, and removed 
before proceeding into attack side decision unit. Except for 


unusual positioning, it isa normally zero for beam type runs. 


4. Iterations begin at Ga) and the time lapse between each suc- 
cessive value of K represents one sampling period. Actual 
program used a scale factor throughout to control the length 


of this period. 
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The following are the subroutines used by the simulator and 
referred to in the block diagram of the system by S/R. 
+ * © oF 
OFFSET subroutine computes Xans and Yo. from entry arguments 


showns 


LOT#2 
Veer = kab, (1- cos(lGel)) 


‘N 


Xeer = Dy 2A0., Sin (| Gp!) 
Site = 169] /Oe 


ats ® Kee, 005 (Org |~ s10€+Yex, Sin (Gre) 


Yabs * Lees Sin (Ore) + $108 *Ypes O05 (2) 





WIND subroutine computes true heading Ty and true airspeed Ty 


by vector subtraction of wind from relative track. ' 
SUBROUTWIE_WINO (Ry LV, We Wy, Te Tv) 












A= Wy -Ry 

B62 Rv+Wv (C05 (6)) 

ied aay (ee Sih e) 
vw 






Tis Rit Pp 
WV = ¥/¢os(p) 













Z 


TURN subroutine finds direction and magnitude of shortest 


turn from P radians to R radians. 


SusRoutwEe TURN (22,6) 
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TRACK subroutine solves ARCTAN problem for correct quadrant 


according to the signs of X and Y arguments. 
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RTRANS subroutine converts radians value G to north or= 


iented direction F in degrees. 


SUBROUTINE RTRANS (G,F) 





TRANS subroutine converts north-oriented heading in degrees 


to radians. 


SUBROUTINE TRANS (D,E) 
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RECTAN subroutine converts polar coordinate to rectangular 


coordinates. 


SUBROUTINE LECTAN (V8. LY) 


Ye Vsin(0) es ae 
; XzV cos (0) 





AERODYN subroutine solves for radius of turn RAD, rate 
of turn ROT, and indicated airspeed for given bank angle @, 


altitude ALT , and true airspeed TAS. 


SuBROUNNE AELOPYN (TAS, G , ALT, RAD, kor, TAs) 







z= tan @ 
RAD = 189.3 (TAS)*/2 
ROT = 005289 Z/ TAS 


IAS = 3600. (FAS)/ (14,02(Atz)) 
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